Competition between pairing and ferromagnetic instabilities in ultracold Fermi gases near 

Feshbach resonances 
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We study the quench dynamics of a two-component ultracold Fermi gas from the weak into the strong inter- 
action regime, where the short time dynamics are governed by the exponential growth rate of unstable collective 
modes. We obtain an effective interaction that takes into account both Pauli blocking and the energy depen- 
dence of the scattering amplitude near a Feshbach resonance. Using this interaction we analyze the competing 
instabilities towards Stoner ferromagnetism and pairing. 



Ferromagnetism in itinerant Fermions is a prime example 
of a strongly interacting system. Most theoretical treatments 
rely on a mean-field Stoner criterion [1], but whether this ar- 
gument applies beyond mean-field remains an open problem. 
It is known that the existence of the Stoner instability is very 
sensitive to the details of band structure and interactions ||2}0J, 
however how to account for these details in realistic systems 
remains poorly understood. Exploring the Stoner instability 
with ultracold atoms has recently attracted considerable atten- 
tion. Following theoretical proposals ||5l, the MIT group made 
use of the tunability (6] and slow time scales ifT UTOl of ultra- 
cold atom systems to study the Stoner instability ifTTIl . Sig- 
natures compatible with ferromagnetism, as understood from 
mean-field theory flTl, were observed in experiments: a max- 
imum in cloud size, a minimum in kinetic energy and a maxi- 
mum in atomic losses at the transition. However, no magnetic 
domains were resolved. 

An impoitant aspect of the MIT experiments is that they 
were done dynamically: the Fermi gas was originally prepared 
with weak interactions and then the interactions were ramped 
to the strongly (repulsive) regime. Dynamic rather than adi- 
abatic preparation was used in order to avoid production of 
molecules. This raises the question of what are the dominant 
instabilities of the Fermi gas in the vicinity of a Feshbach res- 
onance. 

Naively, one would expect that on the BEC-side, molecule 
production is slow, as it requires a three-body process. There- 
fore, instability towards Stoner ferromagnetism would dom- 
inate over the instability toward molecule production. In 
this picture, one would expect that quenches to the attrac- 
tive (BCS) regime always yield an instability towards pairing, 
whereas quenches to the repulsive (BEC) regime an instability 
towards ferromagnetism for sufficiently strong interactions. 

In this Letter, we argue that this picture, which was used to 
interpret the MIT expeiiments, is incomplete. Near the Fes- 
hbach resonance, even on the BEC side, pair production re- 
mains a fast two-body process as long as the Fermi sea can 
absorb the molecular binding energy. As a result, near the 
Feshbach resonance, both on the BEC and the BCS side, the 
pairing and the Stoner instabilities compete directly. We now 
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FIG. 1: Growth rate of the pairing and Stoner ferromagnetic insta- 
bilities after a quench as a function of the final interaction strength 
1/kpa. Final interactions with negative (positive) values of 1/kpa 
correspond to the BCS (BEC) side of the Feshbach resonance. The 
Stoner instability simultaneously occurs in multiple channels. The 
most unstable channel is indicated by the solid red line, the others 
by dashed red lines. The "RPA Stoner" instability corresponds to the 
RPA result with bare as opposed to Cooperon-mediated interaction 
(see text and Ref. fTsl ). 

Inset: Schematic diagram of the pair creation process showing the 
binding energy (spring) being absorbed by the Fermi sea (arrows). 



discuss these instabilities and their competition in detail. 

We start by describing the inter-atomic interactions. A Fes- 
hbach resonance enables tunable interactions between ultra- 
cold atoms by coupling the collision partners to a molecular 
state with different magnetic moment. For broad resonances, 
where the coupling is much larger than the Fermi energy, this 
can be modelled by a single collision channel that supports 
one shallow bound state lfT4l . An often used, but patholog- 
ical choice, is to describe repulsive interactions with a hard- 
sphere pseudo-potential. Instead, one should use the full T- 
matrix that includes the molecular bound state ITSl . Although 
at low energies the scattering amplitude from the hard-sphere 
potential and the T-matiix match, at higher energies compara- 
ble to the molecular binding energy, they do not. Specifically, 
in the strong interaction regime where the Stoner instability 
is expected to occur, the Fermi energy is comparable to the 
binding energy of a molecule in vacuum, precluding the use 
of the hard-sphere potential. 



In light of this remark, we study the initial dynamics of 
the collective modes of a Fermionic system after a sudden 
quench, taking the Cooperon (full T-matrix and Pauli block- 
ing) into account. We focus on the case of a sudden quench, 
as it is simpler and captures the essential physics of the insta- 
bility of the Fermi surface. Extensions to finite rate quenches 
are discussed in Ref. |13|. Our main findings are summa- 
rized in Fig.[T]and are: (1) We find that with the full T-matrix 
the Stoner instability survives with a finite growth rate in the 
range —0.2 < kpa < 1, where a is the scattering length and 
kp is the Fermi momentum. In contrast, using bare interac- 
tions lfT3l results in an unphysical divergence of the growth 
rate at unitarity and no magnetic instabilities on the BCS side 
(see Fig. [T]i. (2) The pairing instability persists on the BBC 
side, where it competes with the Stoner instability. (3) Within 
our approximations, the pairing instability is stronger. 

At first sight, the survival of the pairing and Stoner instabil- 
ities on the wrong side of the resonance is quite remarkable. 
However, both can be understood by taking into account the 
presence of the Fermi sea. On the EEC side, due to Pauli 
blocking, the binding energy of the pair-like molecule can be 
absorbed by the two holes that are left behind (see the inset 
of Fig. [T}. Thus, the two-body pairing process becomes for- 
bidden when the binding energy ^ 1/ma^ exceeds the maxi- 
mum energy that can be absorbed by the holes ^ kp/m (m is 
the Fermion mass, a the scattering length, kp the Fermi mo- 
mentum, and throughout this Letter we use the units in which 
h = 1). On the BCS side, although interactions at low ener- 
gies are indeed attractive, the same is not true at high energies. 
As the Stoner instability involves all scattering energies up to 
the Fermi energy, it is natural that it can persist around unitar- 
ity, even on the BCS side. 

Formalism - We consider a system of interacting Fermions 
described by the Hamiltonian: 



l/k^a 



H^J2 a-4.Ck. + / rf'r U{t, r - r')cl^cl,^ 
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where cl{c^) are the Fermion creation (annihilation) opera- 
tors with spin cr, ^ko- = k'^ /2m — Ha, /Xo- are the chemical 
potentials, and U{t,r — r') is the time dependent pseudo- 
potential that describes the inter-atomic interaction. We fo- 
cus on the instantaneous quench limit, in which the coupling 
U changes from a negligible initial value Ui to a final value 
Uf at time t — 0. In this hmit, we can describe short time 
dynamics of a collective mode at momentum q using the cor- 
responding susceptibility, XqCi^q; Uf), evaluated with final in- 
teractions but initial Fermionic configuration lfT3l[T6l . In par- 
ticular, if Xq(^q; Uf) has poles at Wq = fiq + lAq in the 
upper half of the complex plane, then fluctuations that occur 
after the quench will grow exponentially in time. Next, we 
obtain a universal description of interactions created by the 
pseudo-potential U{r) by modifying the T-matrix formalism 
to take into account Pauli blocking, and apply these ideas to 
the Stoner and BCS instabilities. 

Cooperon - In this section, we obtain the Cooperon, C, i.e. 
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FIG. 2: Pairing instability, (a) "Binding energy" of a Feshbach 
molecule in vacuum and in the presence of a Fermi sea (relative 
to 2eF) as a function of interaction strength, corresponding to the 
real part of the T-matrix pole frequency Qq=o ~ Re[Ljq=o]. Pauli 
blocking by the Fermi sea results in stronger binding across the res- 
onance. The kink occurs when the pair becomes stable, (b) Pairing 
rate and (c) rate of change of Kinetic Energy as a function of interac- 
tion strength on the BEC side for various temperatures [T = (pur- 
ple, solid), 0.12, 0.22, 0.5, 0.66, 0.75 Tf (red, dashed)]. Tempera- 
ture is more effective at suppressing pair production at larger values 
of kfa as the binding energy is smaller, thus the peaks in (b) and (c) 
become sharper at higher temperatures. The peaks in growth rate and 
kinetic energy rate qualitatively match experiments 1 11 1. Sharp onset 
at fc^a ~ 1.1 is expected to be smoothed by three-body processes. 



the T-matrix that takes into account Pauli blocking of states 
by the Fermi sea (see Fig.[3]l, which is needed for an accurate 
study of instabilities near a Feshbach resonance. 

In the center of mass frame in vacuum, the scattering of a 
pair of particles with identical masses m near a wide Feshbach 
resonance is described by the T-matrix (scattering amplitude) 



-(£;) = 



47r 



(2) 



Here, E is the energy of the scattered particles and the pseudo- 
potential U{r — r') that appears in Eq. [T] is related to the 
T-matrix via the Lippmann-Schwinger equation. To cor- 
rectly renormalize the Cooperon, we compare the Lippmann- 
Schwinger equations in a Fermi sea and in vacuum. For the 
Cooperon we can not just use the center of mass frame, as the 
Fermi sea breaks translational invariance. Therefore, we use 
the laboratory frame for both to obtain 



C-\E,<i)^T-^{E + 2e/ - q2 /4m) 
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(3) 



Here, E and q are the center of mass frequency and momen- 
tum of the pair, ef is the Fermi energy, n^ (k) is the Fermi 
function, and et = k^/2m — ef. Our approach is analogous 



to the one used for the Fermi-polaron problem ifTTl . 

Pairing instability - The Cooperon enters the pairing sus- 
ceptibility via 



Xp^M = J dh dh Gih)Giq-h)Ciq)x 

X G{k2)G{q-k2), (4) 

where q stands for the external frequency and momentum 
vector {E,q}, dki stands for doji dki/{2TT)'^, and G(fci) = 
G{bJi, ki) is the bare Fermionic Green function in the non- 
interacting Fermi sea corresponding to the initial state. The 
poles of the pairing susceptibility correspond thus to the poles 
of the Cooperon, whose structure we now explain. 

We begin our analysis with the T-matrix in vacuum. For 
each value of the scattering length, the T-matrix has a line of 
poles on the BEC side located at Wq = i^g+iAg = —l/ma^ + 
mq^/4, corresponding to the binding energy of a Feshbach 
molecule with center of mass momentum q. As a consequence 
of energy and momentum conservation the pole frequency is 
strictly real (Ag = 0), indicating that a two-body process in 
vacuum cannot produce a Feshbach molecule. 

In the presence of a Fermi sea, the states below the Fermi 
surface are Pauli-blocked, shifting the poles of the Cooperon 
relative to the T-matrix in vacuum in two important ways. 
First, the real part of the pole ilg, which would correspond 
to the binding energy of a pair in the absence of an imaginary 
part, uniformly shifts down (see Fig.|2^). This shift is a result 
of Pauli blocking |18|, and indicates an appearance of a paired 
state on the BCS side as well as stronger binding of the pair on 
the EEC side. Second, in the range — oo < l/kpa < 1.1 the 
pole acquires a positive imaginary part Ag that corresponds to 
the growth rate of the pairing instability. As depicted in Fig.[T] 
the growth rate of the pairing instability increases exponen- 
tially Ag=o ~ SeFC^/^''^"^^ as one approaches the Fesh- 
bach resonance from the BCS side, i.e. the growth rate of the 
BCS pairing is equal to the BCS gap at equilibrium 1 18|. On 
the BEC side, the growth rate continues to increase, reach- 
ing a maximum at fc^-a w 2, and finally decreasing to zero at 
kpa w 1.1, at which point the Fermi sea can no longer absorb 
the energy of the Feshbach molecule in a two-body process. 
Pairing deeper in the BEC regime takes place via the more 
conventional three-body process and would round the pairing 
instability curve near kpa ^ 1.1 in Fig.[T| 

So far, we have concentrated on the instability at q ^ 0. 
We find that g = is indeed the most unstable wavevector 
throughout the Feshbach resonance, but the growth rate re- 
mains finite up to q = q^m- Throughout the resonance the 
approximation Qcm ~ i\/S/2){Aq^Q/ep)kp works reason- 
ably well except in the vicinity of /c/a ~ 2 where q^ut reaches 
the maximal value for a two-body process of 2fc/. 

Stone r instability - One can expect that a rapid quench to 
the BEC side of the resonance, where interactions are strongly 
repulsive, results in an instability towards Stoner ferromag- 
netism. We shall assume that right after the quench, the atoms 
are still in the free Fermi sea initial state and the Stoner insta- 




FIG. 3: Diagrams for the vertex function T and Cooperon C. Solid 
lines represent bare fermion propagators, dashed lines interactions, 
and wavy lines external sources. External legs, represented by gray 
lines, are shown for clarity. 




1-0.5 



0.25 



q /k„ 

^max F 

A /^^ 

max F 



0.2 0.4 0.6 0.8 1 10-' 
q/kp 



10-- 



iQ-2 *max F 

10-t A /£„ 

^^ max F 

10-0 



FIG. 4: Properties of growing collective modes in the Stoner insta- 
bility in 3D. (a) Growth rate of the most unstable mode Ag as a 
function of wavevector q for T — and l/kpa — 0.85 (top line), 
0.86, 0.87, ..., 0.93 (bottom line), (b) The most unstable wavevec- 
tor Qniax (blue) and the corresponding growth rate Amax (red) vs. 
l/kpa. A fit to the mean-field critical theory (u = 1/2, z = 3) 
is shown with black lines |19|. (c) Details of the critical behavior 
of gmax and Amax as a function of distance from the transition point 
u = {l/kFa)c - (l/kpa), {l/kpa)c ^ 0.94. 



bility is competing with the pairing instability. Our goal is to 
compute the ferromagnetic susceptibility using the Cooperon 
to describe effective inter-atomic interactions. Using the full 
Cooperon allows us to include three important aspects of the 
problem: energy dependence of the scattering amplitude near 
the Feshbach resonance; Pauli blocking, which renormalizes 
the energy of the virtual two particle bound states involved in 
scattering; and Kanamori-like many -body screening |4|. 

Technically, we compute the vertex function r^^ q(wi, ki), 
which is related to the susceptibility via 



XFM{q)= dkiG{q + k^)G{ki)T^{k^). 



(5) 



We note, that the poles of the susceptibility and the vertex 
function coincide. Replacing the point contact interaction ver- 
tex by the Cooperon in an RPA type resummation of the vertex 
function (see Fig.l3]and Ref. |4|) we obtain 

rg-(fci) = 1 

+ / dk2 T^{k2)C{k +k2+ q)G{k2 + q)G{k2). (6) 



To compute the vertex function, a number of approxima- 
tions are unavoidable. First, we assume that q and ui are 
both small, which is valid in the vicinity of the Stoner tran- 
sition. Second, in the spirit of Fermi liquid theory, we as- 
sume that the most important poles come from the Green func- 
tions, and thus we replace G(k2 + q, a;2 + tj)G(k2,w) — > 
T^TiT^^Ji^)^i\^2\ ~ kp) CD. We then obtain 



where 



Iq_.uj{^2 



kldk2 nF(k2 - q/2) ~ n^ (kz + q/2) 

27r2 UJ - ek2_q/2 + ek2+q/2 



, (8) 



and k indicates a vector on the Fermi surface. We thus make 
the approximation that we can replace ki and kz by ki and kz 
when we evaluate the value of the Cooperon. In other words, 
we assume that the Cooperon changes slowly compared to 
the Green functions. The approximation is fully justified for 
weak interactions, where the Cooperon is momentum and fre- 
quency independent, and the vertex function matches the RPA 
result [ 13 1. For strong interactions, the Stoner instability is not 
driven by the pole of the Cooperon, and we therefore believe 
that our approximation captures the essential physics. 

In the range —0.2 < l/kpa < 1.0, there is one or more 
lines of complex poles with a positive imaginary part Aq, 
which corresponds to the Stoner instability in different chan- 
nels (a combination of momentum and orbital moment). As 
g — > 0, the different instabilities can be identified as different 
angular momentum channels. Since magnetization is a con- 
served order parameter, in each channel Aq grows linearly 
for small q. At large q the cost of bending the order parame- 
ter results in the vanishing of Aq for q > qcut- In between, 
Aq reaches its maximum value Amax at a wave-vector g,nax 
which corresponds to the fastest growing mode (see Fig.Hli. 

Discussion - The growth rates of the pairing instability 
A^£q and the ferromagnetic instabilities in the various chan- 
nels A™^ are compared across the Feshbach resonance in 
Fig. [T] We have also included the naive RPA estimate of the 
growth rate of the Stoner instability in which we have replaced 
the Cooperon by ATra/m. From the comparison, we see that 
(1) the Cooperon suppresses the growth rate of the ferromag- 
netic instability but does not eliminate it, (2) the pairing and 
ferromagnetic instabilities compete over a wide range of in- 
teraction strength on both sides of the resonance, and (3) the 
pairing instability is always dominant. Our results suggest that 
even if there is a metastable ferromagnetic state fTsl, it prob- 
ably cannot be reached via dynamic tuning of the interaction 
starting from a balanced gas. However, for short timescales 
^ (A™^) ^ (A^£q) , both pairing and magnetic corre- 
lations will develop and may be detectable experimentally. 

Comparison with experiment - The maximum of the pair- 
ing instability in the vicinity of fcj^a w 2 closely matches the 
location of the transition found experimentally 1.11.1 . To inves- 



tigate further, we plot the pairing rate as a function of inter- 
action strength for several different temperatures in Fig. Wn. 
The shape of the pairing rate curve, especially at higher tem- 
peratures looks qualitatively similar to the atom loss rate (into 
pairs) found experimentally. 

A fast rampdown of the magnetic field was used to convert 
weakly bound molecules into strongly bound molecules, and 
the kinetic energy of the remaining atoms was measured. It 
was found to have a minimum at fc^a ~ 2 ifTTl . We show that 
this minimum can be qualitatively understood within our anal- 
ysis of the pairing instability. The energy of each molecule 
produced is given by ^ Re[wg] (see Fig.l2k). The molecular 
energy corresponds to the kinetic energy of the Fermions re- 
moved from the Fermi sea, measured with respect to the Fermi 
energy. Thus the rate of kinetic energy change of "unpaired" 
atoms is ~ (Re[cjq] — 2ep) x Im[cjq=o] (see Fig. III). We 
find that the kinetic energy minimum is in the vicinity of the 
maximum of the pairing rate, in qualitative agreement with 
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